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A  STRESS  FACTOR  METHOD  FOR 
PERFORATED  MUZZLE  BRAKE  DESIGN 


Carry  C.  Carofano  and  Martin  R.  Laach 
U.S.  Army  Armament  Research,  Development,  and  Engineering  Center 
Close  Combat  Armaments  Center 
Benet  Laboratories 
Watervllet,  NY 


ABSTRACT 

A  perforated  muzzle  brake  consists  of  a  set  of  vents 
drilled  through  the  wall  of  a  cannon  near  the  musle  to 
reduce  the  impulse  transmitted  to  the  mount.  The  vented 
structure  must  support  both  the  tube  pressure  and  the 
pressure  acting  on  the  vent  surfaces  that  produces  the 
braking  force.  This  report  presents  a  method  for 
estunatmg  the  maximum  Von  Mises  stress  within  the 
structure.  The  results  agree  to  within  10  percent  of  those 
from  a  more  detailed  finite  element  calculation. 


INTRODUCTION 

Perforated  muzzle  brakes  reduce  the  axial  thrust  of  the 
propellant  gas  flowing  from  the  muzzle  by  venting  some  of 
the  gas  through  the  cannon  wall  (see  Figure  1).  The 
pressure  distribution  in  each  vent  is  asymmetric  with  the 
highest  pressure  acting  on  the  downstream  surface,  hence, 
a  net  axial  force  is  developed.  The  left  end  of  the  vented 
region,  called  the  brake  entrance,  must  support  the  total 
brake  load  and  the  maximum  internal  pressure.  In 
practice,  an  iterative  process  is  used  to  And  a  vent  pattern 
which  will  produce  sufficient  braking  without  exceeding  a 
specified  stress  level.  Since  the  flow  and  stress  fields  are 
three-dimensional,  some  simpliflcation  is  needed  to  make 
the  process  practical.  An  approximate  method  for  solving 
the  stress  problem  is  presented  here  that  parallels  the 
approach  used  in  a  previous  analysis  of  the  flow  problem. 

GASDYNAMICS  ANALYSIS 

Earlier  work  (Nagamatsu  et  al.  1987;  Carofano  1988) 
considered  the  three-dimensional  flow  at  the  brake 
entrance.  Because  the  vents  are  placed  symmetrically 


around  the  tube,  the  analysis  was  restricted  to  a  single  vent 
and  the  portion  of  the  tube  volume  associated  with  it.  The 
characteristic  time  for  tube  blowdown  is  determined  by  the 
cannon  length,  while  the  time  to  establish  the  vent  flow  is 
proportional  to  the  tube  wall  thickness,  a  much  smaller  value. 
Therefore,  the  vent  flow  was  treated  as  quasi-steady.  The 
flowfleld  and  pressure  distribution  presented  in  this  section 
were  obtained  using  the  computer  model  described  by 
Nagamatsu  et  al.  (1987)  and  Carofano  (1988). 

The  velocity  vector  and  pressure  contour  plots  in  Figure  2 
show  a  typical  flow  pattern  in  the  sytiunetry  plane  of  the 
vent.  A  uniform  flow  enters  the  brake  at  Mach  number, 
M=  1.5.  and  accelerates  as  a  portion  of  the  gas  expands  and 
turns  into  the  vent.  The  shock  at  the  downstream  lip  of  the 
vent  turns  the  expanded  flow  parallel  to  the  solid  surfaces 
and  reduces  the  velocities  to  subsonic  levels.  The  solid  lines 
in  the  vector  plot  indicate  where  the  local  Mach  number  is 
unity.  The  flow  in  the  tube  accelerates  away  from  this  region 
and  leaves  the  grid  at  supersonic  velocities.  Because  the  tube 
pressure  is  well  above  ambient,  the  flow  is  also  supersonic 
over  most  of  the  vent  exit  area  and  was  found  to  be 
insensitive  to  the  extrapolated  boundary  condition  used  over 
the  subsonic  portion.  For  the  Euler  equations,  then,  the  flow 
is  independent  of  events  occurring  farther  downstream  or  in 
the  surroundings  and  is  completely  described  by  the  vent 
geometry  and  the  Mach  number  and  thermodynamic 
properties  of  the  gas  upstream  of  it.  One  solution  with  these 
parameters  specified  is  valid  for  all  pressures  and  densities. 
Average  values  of  the  mass  flux  at  the  vent  exit  were 
obtained  as  functions  of  the  parameters.  This  information 
was  used  to  construct  a  source  term  in  the  one-dimensional 
time-dependent  Euler  equations  which  relates  the  local  vent 
flow  rate  at  a  specific  axial  location  to  the  local  vent  area  and 
the  pressure  and  density.  The  equations  were  numerically 
integrated  over  the  length  of  the  cannon,  starting  from  a  time 
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when  the  projectile  base  is  just  upstream  of  the  brake 
entrance  and  ending  with  the  completion  of  the  blowdown 
process. 

The  vent  pressure  distribution  is  needed  for  the  stress 
analysis.  It  is  a  function  of  the  vent  geometry  and  the 
local  pressure  and  Mach  number  in  the  tube.  The 
dependence  on  Mach  number  is  shown  in  Figure  3a  tor  a 
geometry  having  a  height-to-diameter  ratio,  h/d,  of  unity. 
The  abscissa,  z/d.  is  the  dimensionless  distance  throu^ 
the  tube  wall  measured  from  the  bore  surface.  The 
ordinate,  p/p.,  is  the  local  pressure  divided  by  the  uniform 
pressure  at  the  brake  entrance.  The  two  curves  in  each 
plot  correspond  to  brake  entrance  Mach  numbers  of  1.0 
and  1.5.  The  higher  Mach  number  produces  a  greater 
pressure  on  the  downstream  surface  behind  the  shock  but 
IS  of  little  significance  elsewhere. 

The  effect  of  vent  geometry,  h/d.  is  confined  to  the 
upstream  wall,  as  shown  in  Figure  3b.  The  pressure  level 
is  higher  for  the  taller  vent  (h/d3:1.25)  because  the  gas 
actually  swirls  around  the  vent  surface  as  it  dows  radi^ly 
outward  and  is  compressed  at  the  upstream  wall  (ref  1). 
Note  the  stagnation  point  in  Figure  2  and  the  recirculation 
zone  below  it.  A  vent  with  h/d=  1  represents  the  crossover 
point  where  the  net  axial  load  on  the  vent  surface  starts  to 
decrease  with  increasing  height. 

There  are  a  nurtiber  of  other  parameters  that  affect  the 
stress  field  and.  in  view  of  the  previous  discussion,  adding 
the  .Mach  number  and  vent  geometry  to  the  list  seems 
unnecessary.  Instead,  the  pressure  distribution 
corresponding  to  M=1.5  and  h/d=l  was  chosen.  Based 
on  internal  ballistics  solutions,  this  Mach  number 
represents  an  upper  limit  at  the  brake  entrance  and  will 
provide  a  conservative  estimate  of  the  maximum  stress. 
Furthermore,  since  h/d=l  represents  the  optimum  dow 
geometry,  it  will  appear  most  often  in  design.  The 
complete  vent  pressure  distnbution  is  tabulated  in 
Appendix  A. 

The  tube  surface  pressure  distribution  near  the  vent  is 
also  needed.  It  is  uniform  upstream  because  the  brake 
entrance  Mach  number  is  supersonic.  Downstream  of  the 
vent.  It  is  nearly  uniform  (see  Figure  2)  and  approximately 
equal  to  the  entrance  value.  To  simplify  the  analysis,  the 
pressure  was  taken  to  be  uniform  over  the  entire  surface 
area  associated  with  the  vent.  Together  with  the  vent 
pressure  distnbution.  this  completely  defines  the 
instantaneous  loads  on  the  brake  for  the  stress  field 
calculation. 


STRESS  ANALYSIS 

Unlike  the  flow  problem,  the  stress  field  near  the  brake 
entrance  depends  on  conditions  downstream  because  the 
total  axial  load  is  transmitted  through  this  region.  /\lso. 
the  abrupt  transition  to  a  perforated  structure  dictates  that 
a  portion  of  the  tube  upstream  of  the  entrance  must  be 
included  in  the  analysis.  While  the  symmetrical  placement 
of  the  vents  around  the  tube  circumference  means  that 


only  one-half  of  one  column  has  to  be  modeled  (columns  run 
axially),  the  computer  memory  for  a  practical  brake  is  too 
large.  As  a  compromise,  two  configurations  were  studied. 
The  fust  is  a  column  model  that  uses  only  a  few  vent  tows  at 
the  brake  entrance  (tows  run  circumferentially)  and  an 
approximate  boundary  condition  at  the  downstream  end. 
This  allows  the  transition  effect  to  be  modeled.  The  second 
considers  a  single  vent  in  a  pattern  of  geometrically  similar 
vents  in  an  infinite  tube.  This  allows  the  effect  of  vent 
diameter  and  spacing  to  be  calculated  economically.  The 
results  are  compared  to  column  model  predictions  to 
determine  if  conservative  estimates  were  obtained.  The  finite 
element  structural  analysis  code  ABAQUS  was  used  for  both 
models  with  a  modulus  of  elasticity  of  206.8  Gpa  and 
Poisson's  ratio  equal  to  0.29. 

The  column  model  is  shown  in  Figure  4a.  Approximately 
14.000  second-order  tetrahedron  elements  (ABAQUS 
C3D10)  were  used,  each  having  10  nodes.  The  left  boundary 
was  located  one  outside  tube  radius  upstream  of  the  brake 
entrance  and  was  constrained  to  move  in  a  fixed  plane 
perpendicular  to  the  tube  axis.  The  right  end  was 
constrained  to  remain  parallel  to  the  left  but  free  to  move 
axially.  Six  vent  rows  were  found  to  be  sufficient  to  limit  the 
influence  of  this  artificial  restraint  on  the  results  near  the 
brake  entrance.  The  surface  pressure  distributions  used  in 
the  calculations  are  specified  later. 

The  single  vent  model  is  shown  in  Figure  4b. 
Approximately  3500  second-order  triangular  prism  elements 
(ABAQUS  C3D15V)  were  used,  each  having  18  nodes.  The 
vent  was  treated  as  one  of  an  infinite  set  equally  spaced  on 
a  tube.  The  left  and  right  boundaries  were  constrained  to 
remain  geometrically  similar  under  load.  Solutions  for  the 
bore  surface  pressure,  the  axial  load,  and  the  vent  pressure 
were  obtained  in  separate  steps.  The  results  are  then 
superimposed  to  obtain  any  desired  loading. 

A  uniform  bore  surface  pressure  was  applied  in  step  I. 
The  boundaries  were  constrained  to  be  planar  and  normal  to 
the  tube  axis.  Equilibrium  was  maintained  by  restricting  the 
left  boundary  to  remain  fixed.  The  nght  boundary  was  free 
to  move  axially. 

In  step  II.  the  boundaries  were  treated  as  above  except 
that  the  axial  displacement  of  the  right  boundary  was 
determined  by  setting  the  integrated  axial  load  on  this  plane 
to  a  specified  value.  In  a  design  problem,  this  is  the  total 
load  generated  by  the  vents  farther  downstream.  The 
instantaneous  load  generated  by  each  vent  is  given  by  the 
one-dimensional  gasdynamics  calculation  of  the  blowdown 
process  (Nagamatsu  et  al.  1987:  Carofano  1988).  This  is 
demonstrated  later  with  an  example. 

The  vent  pressure  was  applied  in  step  III.  Both 
boundaries  were  allowed  to  distort  but  in  such  a  manner  that 
they  remained  geometrically  similar,  consistent  with  the 
infinite  geometry.  To  obtain  a  unique  solution,  additional 
information  was  required  at  the  boundaries.  It  was  desired 
to  have  the  structure  distort  in  response  to  the  vent  pressure 
distnbution  alone.  The  loads  and  displacements  at  the 
boundaries  were  to  be  pan  of  the  solution,  not  part  of  the 
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boundary  conditions.  This  was  achieved  setting  the  sum 
of  the  axial  displacements  of  the  nodes  on  the  left 
boundary  equal  to  zero.  Translatory  rigid  body  motion 
was  also  avoided  with  this  condition. 

Solutions  were  obtained  for  sets  of  the  parameters  listed 
in  Table  1.  With  reference  to  Figure  5.  r^  and  r^  arc  the 
inner  and  outer  tube  radii,  respectively,  d  is  the  vent 
diameter,  and  and  s^  are  the  axial  and  circumferential 
vent  spacings.  respectively,  s,  is  defined  as  the  arc  length 
at  the  mid-wall  radius  (r„-fr^)/2  that  subtends  the  angle 
9—2tr/r\c  associated  with  each  vent,  or 
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Current  design  practice  led  to  the  choices  for  the 
number  of  vent  columns  and  the  wall  ratio.  The  vent 
spacings  were  chosen  because  they  result  in  height-to- 
diameter  ratios  that  span  the  optimum  gasdynamic  value 
of  unitv. 


RESULTS  AND  DISCUSSION 

The  radial,  tangential,  and  axial  stress  components  are 
given  at  the  three  points  on  the  edge  where  the  bore  and 
vent  surfaces  intersect,  designated  A.  B.  and  C  in  Figure 
da.  The  calculated  shear  stresses  at  these  points  were 
found  to  be  negligible.  The  components  are  used  to 
calculate  the  Von  Mises  stress  for  arbitrary  values  of  gas 


pressure  and  axial  loading.  Experience  with  the  column 
model  showed  that  the  maximum  Von  Mises  stress  always 
coincided  with  one  of  these  points.  The  results  are  presented 
in  the  form  of  stress  factors  in  Table  2.  The  negative 
numerical  values,  when  multiplied  by  the  gas  pressure,  p.. 
give  the  compressive  stresses  acting  on  the  tube  and  vent 
surfaces.  The  radial  component  corresponds  to  the  uniform 
pressure  at  the  bore  surface.  The  tangential  and  axial 
components  are  taken  from  the  vent  pressure  distribution 
tabulated  in  Appendix  A. 

The  factor  is  defined  as 


The  stress.  <t.  is  the  sum  of  the  contributions  from  steps  I 
and  III  that  result  from  the  gas  pressure  loading.  The 
denominator  is  the  tangential  stress  at  the  bore  surface  of  an 
unvented  tube. 

The  factor  iC(^  in  Table  2  is  associated  with  the  axial 
loading  in  step  II  and  is  defined  as 


The  denominator  is  the  axial  stress  in  an  unvented  tube 
having  the  load.  L  distributed  uniformly  at  each  end. 


Table  1.  Dimensionless  Parameters  Describing  Vent  Geometry 


number  of  vent  columns, 

10  and  12 

wall  ratio.  r,T^ 

1.35  and  1.45 

axial  spacing,  d/x^ 

0.50.  0.67  and  0.80 

circumferential  spacing,  d/s^ 

0.47,  0.55  and  0.66 

Table  2.  Stress  Factors 


Location 


Radial 


Gas  Pressure  Loading 


Tangential 


Axial  Loading 

Radial 

Tangential 

Axial 

0.0 

K,^ 

0.0 

0.0 

0.0 

Klb 

0.0 

0.0 

¥if  and  K(.  are  given  in  Figure  6  for  a  geometry  having 
twelve  vent  columns  and  a  wall  ratio,  r,/r^sl.3S.  The 
abscissa  in  each  plot  is  the  circumferential  spacing 
parameter,  d/s^  Each  curve  is  for  one  value  of  the  axial 
spacing  parameter,  The  vent  section  in  Figure  6c 
shows  that  stress  factor  applies  at  point  A,  which  is 
located  on  the  downstream  lip  of  the  vent  on  the  bore 
surface,  and  similarly  for  the  other  factors.  Because  the 
loading  is  symmetrical  in  step  II.  K^;  is  equal  to  and 
is  not  plotted. 

The  trend  of  the  data  is  best  understood  by  considering 
the  load  bearing  area  or  web  between  vents  relative  to  the 
unvented  tube  (see  Figure  4a).  For  the  circumferential 
direction,  this  is  the  ratio  of  the  shaded  area  to  total  area 
in  Figure  5.  or 


cinumferenttai  web  ana  _  ^  fl-— ] 

total  area 

(4) 


As  d/x,  increases,  the  web  becomes  thinner.  Since  the  gas 
pressure  and  axial  load  were  fixed  in  the  analysis,  the 
tangential  stress  factors  Kk.  and  tise  as  d/x^ 
increases  in  magnitude  from  the  lower  to  the  upper  curves 
in  each  plot.  The  area  ratio  is  independent  of  d/s^  which 
explains  the  weak  'dependence  on  this  parameter. 

For  the  axial  direction. 


axuU  wet  area 
total  area 


(5) 


where  Eq.  ( 1 )  was  used  to  simplify  the  expression.  The 
web  becomes  thinner  as  d/s^  increases  and  the  axial  stress 
factors  iC^  and  increase  accordingly.  The  area  ratio 
IS  independent  of  d/x^.  so  the  weak  dependence  on  this 
parameter  is  expected. 

The  results  for  the  other  combinations  of  the  wall  ratio, 
r,  r.,,  and  the  number  of  vent  columns,  n^,  do  not  differ 
appreciably  from  those  in  Figure  6  and  are  not  plotted. 
Instead,  they  are  given  in  polynomial  form  in  Appendix  B. 

The  finite  element  solution  for  step  III  is  shown  in 
Figure  7  for  a  geometry  having  a  thin  circumferential  web. 
The  view  is  from  inside  the  tube  looking  at  the  bore 
surface.  The  vent  surface  is  loaded  asymmetrically  by  the 
gas  pressure  in  this  step  and  the  thin  web  undergoes  a 
large  distortion  near  the  bore  surface.  Interestingly,  the 
tangential  stress  at  point  C  exceeds  that  at  point  A, 
although  the  latter  sees  the  high  pressure  behind  the 
shock.  This  can  be  undeistood  by  recalling  that  the  vent 
IS  one  of  an  infinite  set  so  that  C  also  lies  in  the  next  vent 
downstream,  on  the  same  web  with  A.  Both  points  are 
placed  in  tension,  but  the  stress  at  C  is  higher.  This  is 
true  for  all  of  the  thin  web  geometries,  as  shown  by  the 
upper  curves  in  Figure  6a.  The  stress  factors  and  Kfc 
include  contributions  from  steps  I  and  II.  but  these  steps 


produce  symmetrical  results.  The  difference  between  them 
is  due  to  the  asymmetrical  loading  in  step  HI. 


COMPARISON  OF  THE  TWO  MODELS 

The  gas  pressure  decreases  with  axial  distance,  x.  through 
the  brake  due  to  venting.  The  comparison  is  made  for  a 
pressure  of  unity  at  the  brake  entrance  and  a  pressure 
gradient  equal  to  S  percent  per  vent  through  the  vented 
section.  The  column  model  has  six  rows  of  vents  with 
uniform  spacing.  Xc  and  an  unvented  section  at  the  left  of 
length,  r^  so  the  axial  pressure  distribution  is 

p,  »  1.0  ,  0.0  s  X  s  r 

(6) 

p,  =  1.0  -  0.05(x-r^x,  ,  X  z.  r^*6x^ 

The  vent  pressure  distribution  in  row  number  n  is  based  on 
the  average  pressure,  p...  acting  over  the  length  x,  and  is 
given  by 

p„  =  1.0  -  0.025(2«-l)  .  a  =  1,  6 

To  estimate  the  stresses  using  the  single  vent  model,  the 
axial  load.  L«,  applied  to  the  downstream  boundary  of  vent 
row  number  n  is  written  as 

L.  •  (1/2)  E  Is,  .  »  ■  1.  5 

•■•♦I 


At  the  brake  entrance,  n  is  1  and  the  five  vents  farther 
downstream  each  contribute  an  axial  load  L...  At  the  sixth 
vent.  L,  is  zero,  since  no  load  was  applied  at  the  right 
boundary  in  the  column  model.  The  factor  1/2  appears 
because  the  column  model  is  made  up  of  half-vents.  U.,  is 
the  integral  of  the  vent  pressure  distribution  at  vent  number 
m  and  is  shown  in  Appendix  A  to  be 

=  0.9  (rtdV4)  (9) 

p„  is  the  average  pressure  acting  over  the  length  x^  of  vent 
m  and  is  given  by  Eq.  (7). 

Since  the  shear  stresses  are  negligible  at  the  three  points 
of  interest,  the  Von  Mises  stress,  may  be  written  as 


where  and  a,  are  the  radial,  tangential,  and  axial  stress 
components,  respectively.  These  are  given  directly  for  the 
column  model  by  the  finite  element  solution.  For  the  single 
vent  model,  they  are  computed  using  the  stress  factors  in 
Table  1 
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The  following  set  of  parameters  was  used: 
rjr^  »  1.35  ,  •  12 

dls^  =  0.55  .  d/x^  =  0.fi7 


(11) 


From  Figure  6. 


’  2.76 

.  =  -0.05 

=  0.26 

.  =  3.03 

(12) 

Kfc  =  2.75 

,  =  -0.05 

From  Table  2,  the  stresses  for  points  A.  B.  and  C  in  vent 
number  n  are.  respectively. 

=  -P« 


n  ^  JP 

f  1  V  1 

Wjrf-y 

Pm 

13) 


o,  =  -2.208p. 


=  -Pm  ' 

<3,  =  -0.163p^  (14) 


r 

I  *'1, 

*  IT  (  ‘  1 

V 

-Pm 

r 

-ir  f  ^  1 

^"1 

o,  =  -O.OSlp^ 


Eqs.  (1).  (8),  and  (9)  may  be  combined  to  give 


i  *  l^cj 

[  -c 

r  Pm, 


(16) 


Eqs.  (12)  through  (16)  express  the  stress  components  in 
terms  of  the  four  dimensionless  parameters  and  the  local 
vent  pressure. 

The  Von  Mises  stresses  are  plotted  in  Figure  8  for  a 
geometry  having  comparable  web  thicknesses,  as  indicated 
by  the  vent  figure  in  the  upper  right-hand  corner.  The 
solid  and  dashed  lines  represent  the  single  vent  and 
column  models,  respectively. 


For  the  column  model,  the  unvented  section  at  the  left 
provides  support  near  the  brake  entrance  that  significantly 
reduces  the  stresses  in  the  first  vent  at  points  A  and  C.  At 
the  right  end.  the  results  for  vent  6  are  obviously  influenced 
by  the  boundary  approximation  and  should  be  ignored. 

The  maximum  stress  in  the  colurrm  model  occurs  at  point 
A  in  the  second  vent.  Because  of  the  nature  of  the  single 
vent  model,  the  maximum  stress  will  always  appear  in  the 
first  vent.  For  this  geometry  and  loading,  the  maximum  is 
correctly  predicted  to  lie  at  point  A  and  is  conservatively 
estimated  with  respect  to  the  column  model,  exceeding  the 
maximum  in  the  latter  by  10  percent.  By  'conservatively 
estimated,'  it  is  meant  that  the  single  vent  prediction  in  the 
first  vent  exceeds  the  maximum  stress  occurring  anywhere  in 
the  column  model.  This  is  the  Von  Mises  stress  that  would 
be  used  in  a  failure  criterion  to  check  structural  integrity. 

For  both  models,  the  stress  at  point  B  is  always  highest  at 
the  brake  entrance  where  the  total  axial  load  must  be 
supported.  The  maximum  stress  in  the  column  model 
exceeds  the  single  vent  estimate  by  S  percent.  If  an  axial 
load  is  imposed  at  the  downstream  end  of  either  model,  to 
simulate  the  addition  of  mote  vents,  the  stress  at  point  B 
simply  increases  uniformly  along  the  brake.  The  discrepancy 
at  the  brake  entrance  remains  the  same,  although  the  stress 
at  point  B  could  become  the  maximum  stress  in  the 
structure.  For  this  geometry,  tne  error  is  well  within  any 
reasonable  choice  for  a  safety  factor. 

To  determine  if  the  situation  might  worsen  for  one  of  the 
other  geometries,  two  extremes  were  examined.  The  first 
considers  the  configuration  having  the  thinnest  axial  web  and 
the  thickest  circumferential  web.  As  shown  in  Figure  9.  the 
maximum  stress  occurs  at  point  B  and  is  conservatively 
estimated.  At  the  opposite  extreme,  the  thinnest 
circumferential  web  and  the  thickest  axial  web.  the  maximum 
stress  occurs  at  point  C.  as  shown  in  Figure  10.  and  is  again 
conservatively  estimated. 


SUMMARY  AND  CONCLUSIONS 

A  stress  factor  method  was  presented  for  estimating  the 
stress  field  in  a  perforated  muzzle  brake.  The  maximum  Von 
Mises  stress  was  found  to  occur  at  one  of  three  critical 
locations  depending  upon  the  geometry  and  loading 
conditions.  Compared  with  the  results  from  a  mote  detailed 
finite  element  calculation,  the  method  predicts  the  location 
correctly  and  estimates  the  maximum  stress  to  within  10 
percent.  The  iteration  time  to  find  an  acceptable  design  is 
significantly  reduced.  The  final  design  can  be  checked  using 
the  more  detailed  model  to  assure  structural  integrity. 
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APPENDIX  A 

The  vent  pressure  distribution.  p(9.Z)/pv  used  in  the 
analysis  is  given  in  Table  Al.  p(8.Z)  is  the  local  wall 
pressure,  p,  is  the  uniform  pressure  at  the  brake  entrance. 
Zsz/d  is  the  dimensionless  distance  through  the  wall 
measured  from  the  bore  surface,  and  d  is  the  vent 
diameter.  9  is  zero  degrees  on  the  downstream  lip  of  the 
vent  (point  A  in  Figure  6c). 

The  integral  of  the  vent  pressure  distribution.  L..  is 

(Al) 


APPENDIX  B 

The  stress  factors  are  given  the  following  polynomial: 
r  -  3.0  ♦  c,(4rx^  ♦  c.^dixf  *  (bd 

where 

C,  ■  c,  ♦  c^dtsj  * 

C,  -  c*  ♦  e^dlsjl  *  c^dlsf 
’  Cj  *  c^dls^  *  c^dls/ 

and 

e,  »  kj^i  *  *  1.  9 

The  coefficients,  k^.  ms  1.18.  are  given  in  Table  B1  for 
n^s  10  and  Table  B2  for  n^s  12.  As  an  example,  for  a  brake 
having  ten  columns  of  vents,  the  stress  factor  is 
computed  from  Eq.  (Bl)  using  the  k’s  in  the  second  column 
of  Table  Bl. 


Tabhi  Al.  V«nt  Praasur*  Distribution,  p(d,Z)/p.. 


TabI*  B1.  Coaffldwit*,  k..  In  Eq.  (B1)  tor  n.alO. 


9.752 

-5.042 

-59.674 

67.375 

24.631 

-35.992 

-33.002 

13.163 

206.175 

-235.675 

-122.999 

162.179 

20.867 

1.151 

-167.704 

195.380 

117.576 

-156.418 


204.136 

204.160 

-239.545 

-94.471 

149.518 

415.876 

-516.J42 

-459.704 

562.812 

61.053 

-91676 

-256.905 

320.570 

275.544 

-340.425 


20.249 

-22.639 

-49.987 

75.967 

11.010 

-38.702 

-88.450 

92.583 

178.664 

-263.169 

-55.611 

139.403 

85.096 

-84.260 

-161.741 

234.663 

57.362 

-130.055 


-17.097 

0.900 

-23.002 

13.057 

43.786 

-47.052 

44.462 

-19.943 

28.687 

3.032 

-106.586 

107.092 

-36.848 

27.252 

19.808 

-50.472 

42.909 

-33.577 


116.119 

-124.265 

-535.074 

571.614 

588.101 

616.450 

-143.995 

134.020 

744.062 

-735.554 

-889.515 

886.312 

42.308 

-27.959 

-271.337 

235.248 

367.002 

-338.712 


Tabto  B2.  Coafflelanta,  l(„,  In  Eq.  (B1)  tor  n,a12. 
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FIGURE  1,  A  PERFORATED  MUZZLE  BRAKE. 


FIGURE  2.  PRESSURE  CONTOUR  AND  VELOCITY  VECTOR  PLOTS  FOR  A 
SINGLE  VENT  AT  THE  BRAKE  ENTRANCE  THE  SOUD  UNES 
IN  THE  VECTOR  PLOT  INDICATE  A  MACH  NUMBER  OF  UNITY. 
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(a)  Effect  of  Mach  number  on  vent  pressure  for  h/d- 


FIQURE  3.  VENT  PRESSURE  DISTRIBUTIONS 


FIGURE  4  THE  COLUMN  MODEL  (a)  AND  THE  SINGLE  VENT  MODEL  (b). 

THE  BORE  SURFACE  IS  AT  THE  BOTTOM  OF  EACH  MODEL. 


tuba  cantarllne 


FIGURE  5.  VENT  NOMENCLATURE.  THE  WEB  IS  THE  TOTAL 
SHADED  AREA  IN  EACH  SECTION. 
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(a)  Streaa  factors  for  gas  pressure  loading. 
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FIGURE  6.  STRESS  FACTORS  FOR  =  1,35  and  n, 


0.00  5.00  10.00  15.00  20.00  25.00 


FIGURE  7  SINGLE  VENT  MODEL  SHOEING  CIRCUMFERENTIAL 
WEB  DISTORTION.  PARAMETERS  ARE  1  35. 
=  12.  d/a.  *  0.47.  AND  d/x,  »  0.80. 


FIGURE  8.  COMPARISON  OF  THE  VON  MISES  STRESS  FOR  BOTH 
MODELS.  VENT  NUMBER  1  IS  AT  BRAKE  ENTRANCE 
tjx^  =  1.35.  n,  =  12.  d/s,  =  0.55.  and  d/x,  =  0.67, 
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Oy/Pe  OV/P* 

0.00  0.00  10.00  10.00  20.00  20.00  0.00  0.00  10.00  10.00  20.00 


-  Column  modol 

-  Slnglo  vont  a»dol 


- r  I  I  I  I - 1 - 

1  2  3  4  S  6 

v«nt  ifuBbor 

FIGUBE  9.  COMPARISON  OF  THE  VON  MISES  STRESS  FOR  BOTH 
MODELS.  VENT  NUMBER  1  IS  AT  BRAKE  ENTRANCE. 
rj\  =  1.35,  =  12.  d/s^  =  0.66,  and  d/x^  =  0.50. 


FIGURE  10.  COMPARISON  OF  THE  VON  MISES  STRESS  FOR  BOTH 
MODELS.  VENT  NUMBER  1  IS  AT  BRAKE  ENTRANCE. 
=  1.35,  Og  =  12,  d/Sj  =  0.47,  and  d/Xg  =  0.80. 
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